The effect of turbulence on the particle impaction on a cylinder in a cross flow 
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Particle impaction on a cylinder in a cross flow is investigated with the use of Direct Numerical 
Simulations (DNS) and with a focus on the effect of turbulence on the impaction efficiency. The 
turbulence considered is isotropic homogeneous turbulence with varying integral scales. It is found 
<p«. that for particles with Stokes numbers in the boundary stopping mode there is up to 10 times more 

^ front side impaction for turbulence with a large integral scale than for a corresponding laminar flow. 

For decreasing integral scales the effect of the turbulence on front side particle impaction efficiency 
is decreasing. The back side impaction efficiency is also found to be influenced by the turbulence, 
but for the back side the strongest effect, and largest impaction efficiency, is found for small integral 
scales. 

V 

^ I. INTRODUCTION 

5/3 Particle-laden fluid flows are common both in nature and in a large number of industrial applications. Depending 

O on the particle and fluid properties, particles may impact on a solid object entrained in the fluid flow. Such particle 
c/) impaction may either lead to the buildup of a deposition layer on the solid-fluid interface, or to erosion of the solid 

object. Typical industrial applications are filters, furnaces, industrial boilers, fan cooled electronics and ventilation 

systems. 

I For several of the applications mentioned above the solid object may be approximated by a cylinder. Due to this, 

and to the simplicity of the cylindrical geometry, determination of particle impaction efficiency on a cylinder in a cross 
T-H flow has become a benchmark case. By using potential flow theory for the fluid flow particle impaction efficiencies 
can been found [IHS] as a function of the particle size. Since potential flow theory assumes infinite Reynolds numbers, 
but still does not account for turbulence, it is clear that this method has major shortcomings. It is nevertheless 
a well accepted method for Stokes numbers larger than unity. The Stokes number St = ''p/'''/ is the ratio of the 
I particle Stokes time to the timescale of the fluid flow around the solid object. Numerous other approaches focusing 
on experimental [l^'S], numerical CFD [HHS] and phenomenological modelling f9]|T0] are also found in the literature. 

As the Reynolds number is increased turbulence will become important and for very large Reynolds numbers 
turbulence will significantly alter the result. Turbulence is either generated due to strong shear induced by the 
cylinder, or, alternatively, the flow approaching the cylinder is already turbulent. For Re = UqD/v ~ 10^ the wake 
downstream of the cylinder will break up into turbulence even when the flow is laminar upstream of the cylinder. 
For such low Reynolds numbers turbulence is not expected to have any impact on the particle impaction efficiency as 
the transition to turbulence appears too far downstream of the cylinder. Here Uq is the mean flow velocity, D is the 
cylinder diameter and v is the kinematic viscosity. As the Reynolds number increases, the transition to turbulence 
^ moves upstream towards the cylinder and for Re ^ 10^, the boundary layer around the cylinder will be turbulent. 

The effect of turbulence on particle impaction has been studied by Douglas & Ilias (1988) ^H,. In this study the 
cylinder was situated within a channel such that the turbulence was generated by the turbulent channel flow and not 
by the cylinder itself. Except from this, not much is found in the literature on the effect of turbulence on the impaction 
efficiency on a cylinder in a cross flow. There is however a large number of publications on particle impaction on the 
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walls of a turbulent channel flow, where the particles are impacting on the channel walls due to turbophoresis (see 
e.g. Picano et al. |T2] and references therein). For a cylinder in a cross flow, turbophoresis is expected to be relevant 
for very large Reynolds numbers Re > 10^. 



II. THE EFFECT OF TURBULENCE ON IMPACTION EFFICIENCY 



As already mentioned, turbulence can be introduced in two different ways; either the turbulence is generated due 
to the shear induced by the cylinder, or turbulence is generated somewhere upstream of the cylinder and is inherent 
in the flow as it approaches the cylinder. The current work focuses only on turbulence generated upstream of the 
cylinder, and does not consider Reynolds numbers large enough for turbulence to be generated in the boundary layer 
around the cylinder. Furthermore it is assumed that the turbulence approaching the cylinder is homogeneous and 
isotropic. 

If the turbulence is generated upstream of the cylinder there are several possible scenarios, some of which will be 
explained in the following. 



A. Large eddy Stokes numbers 

The eddy turnover time is given as 

Te d/Urnis, (1) 

where d is the integral scale of the turbulence and Urms = \/ {{U — UqY) is the root mean square velocity of the 
turbulence and () symbolizes ensemble averaging. If the eddy turnover time is much shorter than the particle Stokes 
time, i.e. Steddy = Tp/T^ » 1, the particles are too slow to be affected by the turbulent eddies. When this is the 
case the turbulence has no effect on the particles, and consequently the impaction efficiency is not affected. Here the 
particle Stokes time is 

(2) 

where S — Pp/ p, p is the fluid density, pp is the particle material density and dp is the particle diameter. 



B. Large scale turbulence 

It is well known that the impaction efficiency 77 = A^impact /Mnit in a laminar flow is determined by the Stokes 
number |13j . Here the number of particles impacting on the cylinder surface is given by Mmpact while iVinit represents 
the number of particles initially positioned such that the mean flow at their position is moving in the direction of the 
cylinder. When the eddy turnover time is equal to or larger than the particle Stokes time, i.e. Tg > t^, the particles 
will follow the turbulent eddies. Consequently, the particle velocities will deviate from the mean flow velocity. When 
the scale of the turbulence is not small compared to the size of the cylinder this yields a modified Stokes number, 
which will then also give a change in the particle impaction efficiency. 

The Stokes number St = TpjTf is proportional to the mean fluid flow velocity since tj = D/Uq. In the turbulent 
case the magnitude of the fluid flow velocity [/, in general different from the mean flow velocity Uq, is stochastic. 
When the integral scale of the turbulence is large, i.e. d ^ D, U may be seen as the mean flow velocity at that 
instant in time. Thus, the Stokes number also is a stochastic variable, effectively being different from the 'laminar' 
St, expressed by the laminar fluid velocity Uq. The effective Stokes number can be expressed as 

St 

Steff =^U. (3) 



As Steff is a linear function of its variance becomes 
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yar (Steff) = (T^t = ) Var{U), (4) 
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cf. Walpole et al. (2007) [H]. Since Uq is constant, ^ 
the Stokes number equals 



shows that crg^ 



£;(steff) = St, 



when U = Ua- The expectation value of 



(5) 



since U fluctuates symmetrically around Uq. 

With a fluctuation in U, the effective Stokes number becomes Stcft = St + A, with A — StSjj/Uo being the resulting 
fluctuation in the Stokes number when 6u ~ U — Uq. Thus, a Taylor expansion in the small parameter A yields, by 
using ([5]), 



?7(Stcff) 



77(St + A) 
77(St) + r;'(St)A- 



77" (St) 



(6) 



The expectation value of this becomes 



£;[?7(Steff)] = 77(St)+77'(St)^[A 
« 77(St) 



V'(st) 



77"(St) 2 
^^'^Sf 



(7) 



By definition, a^^ = ElA"^] — (i?[A])^ — £'[A^] since E[A] = due to the symmetry of the velocity fluctuations 

around the mean. Furthermore, higher order terms, 0{E[A'^]), have been neglected. When there is no turbulence, 
such that A = 0, the expectation value of the impaction efficiency is -E(77(Stefr)) — 77(St), i.e. it equals the impaction 
efhciency in the laminar expected. 

In order to use Eq. Q to obtain the expected values of the front side impaction efficiencies with turbulence present, 
values for 77" (St) are required. 

Furthermore, cr|^ is needed to determine the expectation value of the front side impaction efficiency ?7(Steff) for a 
given St. For non-inertial particles, i.e. St — >■ 0, particles follow the fluid flow exactly and (t|^ = Wrms- For non-zero 
Stokes numbers the particles do not follow the fluid flow and cr|^ must be found from the simulations. This is done 
by looking at how the velocity of the particles deviate from the mean flow velocity Uq. This deviation expresses how 
the turbulence affects the particles, and will typically be different for different Stokes numbers. To find reliable values 
of CTg^, the position where is measured, should not be too close to the cylinder. 



C. Small scale turbulence 



If the turbulent eddies are very small, they may even penetrate into the boundary layer around the cylinder. If this 
happens, the particles can impact on the cylinder surface due to turbophoresis, which may have a significant effect 
on the impaction efficiency. Let us call this effect impaction by external turbophoresis. The boundary layer thickness 
around the cylinder is given by |13j 

VRe 

where B 0.5 such that in order to have impaction by external turbophoresis the turbulent integral scale d must 
fulfill the inequality 

d<— (9) 
2%/Re 



The dissipative term, in spectral space, is given as vkF'(U — Uq) such that for small scale turbulence, i.e. when the 
wavenumber k is large, the decay will be very fast. Consequently, if the turbulence generated upstream of the cylinder 
shall survive all the way down to the cylinder, the timescale of the eddy decay, given by t^^ = (P /ly, must be longer 
than the typical fluid timescale t/ = D/U such that 



1 /Z)^ ^ 



Re>^ -). (10) 
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By combining Eq. ^ and Eq. ( 10 ) it is found that 

Re>^^) >4Re. (11) 

It is clear that the above inequahty is a contradiction, and consequently the conclusion can be drawn that impaction 
by external turbophoresis, where the scale of the turbulence is small enough to penetrate the boundary layer around 
the cylinder, is not practically feasible if the source of turbulence is not very close to the cylinder. An example where 
the source is indeed close to the cylinder is a cylinder placed in a turbulent channel flow [11] . This is a situation which 
is encountered in many industrial applications, but it does introduce some extra parameters into the study. Due to 
the increased parameter space this application is not considered here; instead, this work focuses solely on isotropic 
decaying turbulence. 



III. SIMULATIONS 



The DNS code used for the simulations performed in this work is The Pencil Code [TS]. The fluid flow is evolved 
in time by solving the continuity equation 



dp 

dt 



V • (up) = (12) 



and the momentum equation 



du 11 F 

— +M • Vm+ -VP = -V • (2pi^S) + — (13) 
ot p p p 

together with the equation of state given by P = pc^. Here u is the fluid velocity, Cs is the speed of sound, P is the 
pressure, F is some external force and S is the traceless rate of strain tensor. 
The particles are tracked individually in a Lagrangian manner as 

and 

when V is the particle velocity, x is the particle position and Fp = {u — v)/Tp is the drag force from the fluid on the 
particles. 

In the following, all variables are non-dimensionalized by the mean flow velocity J7o, and the length of one of the 
sides of the simulation box Pbox- 

Isotropic turbulence is created inside a cubic domain, the turbulent box. The length of each of the sides of the box is 
-^box = 1 and two different Reynolds numbers are simulated, namely 420 and 1685. At Re = 420, the number of grid 
points are 512^, while at Re — 1685, 1024'^ grid points are required. Statistically stationary turbulence is achieved by 
the use of external forcing. Energy is added at a given wave number, and the turbulence develops until the energy 
input through the forcing is equal to the energy dissipated in the dissipative subrange. This leads to an equilibrium 
where fluid motion is independent of the turbulence's initial conditions [Hj and the turbulent flow is statistically 
stationary. The forced turbulence is created by adding a stochastic forcing function 

F{x, t) = 3fi{7Vffc(t) exp[ifc(i) • X + i(t){t)]) (16) 



as the external force on the right hand side of the momentum equation (Eq. (13l). As k{t) and <f>{t) are chosen 
randomly at each time step, the stochastic nature of turbulence is inherent in the equation. The normalization factor 
is A'^ = /oc(|fc|c/At)^/^ and ^k(i) is perpendicular to k and an eigenfunction to the curl operator. A forcing wavenumber 
kf is chosen such that the magnitude of the randomly chosen k{t) is in the range kf — 0.5 < \k{t)\/kQ < fcf + 0.5, where 
ko = 1/Pbox is the normalizing wavenumber corresponding to the side lengths Lbox of the box. For more details on 
the forcing see e.g. Haugen & Brandenburg (2006) [17]. The effect of the forcing is that turbulent energy is put into 
the system at the spatial scale corresponding to the forcing wavenumber fcf, which in this manner determines the 
behaviour of the turbulence. As k{ is normalized by fco, the forcing length scale will he d = L^ox/kf = fc^^. 
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(b) Turbulent inlet on the domain. 



FIG. 1: Figure (a) shows how turbulent velocity information containing strips (yellow and blue here) at discrete 
increasing x-values are successively extracted from the xy-plane, i.e. for constant z. As all strips for z = are 
extracted, in the blue marked plane, the procedure is repeated at the next discrete z-value. When the last xy-plane, 
the orange one, is used, the procedure starts over again. In (b) it is shown how the turbulent velocities of the 
successive strips are imposed on the boundary of the domain. 



The turbulence is imposed on the inlet of the two-dimensional main domain (see Fig. |lb[) of size L^; x = 2 x 1. 
The cylinder cross section with a diameter D = Lbox/6 is placed in the center of the 2D domain. The fluid flow enters 
the domain with a mean velocity Uq ^ 1. Prior to letting the turbulence enter the domain, von Karman vortices in 
a steady state are established in the wake of the cylinder. When this is done, the turbulence, generated in the box 



shown in Fig. la is used as inlet for the two dimensional domain containing the flow, the cylinder and the particles. 
As illustrated in Fig. [TJ a quadratic xy-slice of the turbulent box is divided into strips whose velocity information is 
extracted and inserted at the leftmost side of the domain, shown in Fig. |lb[ At successive time steps, the position 
of the strips chosen depends on the distance traveled by the fluid, i.e. U^t. When the end of a slice is used as inlet, 
the slice at the next discrete z-value is used in the same manner. When the outer end of the box is used, i.e. the 
slice at the maximum z-value, the procedure is started over again with the strips from the slice at z = 0. Imposing 
turbulence in this way essentially means adding the turbulent velocity Ut to the velocity UqX. Hence, the velocity at 
the inlet is 

Mi„ = C/o* + Mf (17) 

The boundary conditions implemented when imposing the turbulent velocity at a; = are the Navier-Stokes charac- 
teristic boundary conditions (NSCBC) [18]. The advantage of using NSCBC is that the boundaries are non-reflecting, 
meaning that any signal is let through them. The direct numerical simulation of compressible flows requires an ac- 
curate control of wave reflections from the domain boundaries, as the accuracy of the solution is in general sensitive 
to solutions at the boundaries. In compressible fluid flows, and especially in DNS where the range of scales is large, 
reflected waves can cause problems. This is due to e.g. convected vorticity or sound waves not being let through 
the boundaries, but rather being reflected back into the domain. This may even lead to strong standing waves in 
the computational domain. The transparency of NSCBC prevents this. In the current work the NSCBC method as 
described by Lodato et al. (2008) [1^ is used both for the inlet and the outlet. 

The 2D domain containing the cylinder cross section and the surrounding flow, where the flow is in the x-direction, 
has periodic boundary conditions in the y-direction, meaning that a particle or a fluid element hitting the boundary 
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at y = or y = Ly immediately is reinserted at the same x-position on the other side, with the same velocity and in 
the same state. 

Particles are inserted at a; = 0, and are removed from the simulation when hitting the cylinder or the rightmost 
boundary at x = L^. Among similar simulations done in the past, there have been some disagreements related 
to the correct number of particles to insert to achieve the desired statistical reliability of the data [5D]. To ensure 
statistically significant data, a large number of particles, N — 10^, is inserted in each simulation. The particles are 
inserted with an initial velocity Vq = J/q at a rate of 0.133 • 10^ particles per time unit, where the time unit has been 
non-dimensionalized hy t = Li^ox/Uq. This correspond to a total particle insertion time of several von Karman eddy 
times, which is required in order to get statistically stationary results. An overview of the two-dimensional simulations 
is given in table |l] 

TABLE I: Overview of the 2D simulations. 'L' corresponds to lower Re runs and 'H' to higher. 



Case 


Resolution 


Re 


Flow regime 


LI 


1024 X 512 


420 


Laminar 


L2 


1024 X 512 


420 


Turbulent; fcf = 15 


L3 


1024 X 512 


420 


Turbulent; fcf = 5 


L4 


1024 X 512 


420 


Turbulent; fcf = 1.5 


HI 


2048 X 1024 


1685 


Laminar 


H2 


2048 X 1024 


1685 


Turbulent; fcf = 15 


H3 


2048 X 1024 


1685 


Turbulent; k{ — 5 


H4 


2048 X 1024 


1685 


Turbulent; k{ = 1.5 



IV. RESULTS 



The relative difference in impaction efficiency between a turbulent and a laminar simulation is defined as 

R=R^3^ (18) 

where r] is the turbulent impaction efficiency and 771am is the impaction efRciency of the corresponding non-turbulcnt 
setup. In Fig. [2] the relative impaction efficiency is shown for simulations with Re = 420 and for forcing aX kf — 1.5 
and kf — 5 both for two and three dimensional simulations. It is clear that there are some differences between the two 
and the three dimensional results, but qualitatively the results are similar. In the following the focus will therefore 
be on the two dimensional results as three dimensional simulations on Reynolds numbers larger than 420 can not be 
afforded. 



A. Effect of turbulent integral scale for Re — 1685 



Figure |3] shows the front side impaction efficiency plots, with their characteristic shape, for simulations with Re — 
1685. The boundary stopping mode is in the range 0.15 ^ St < 0.5, with the inertial impaction and boundary 
interception modes above and below this, respectively. The effects of turbulence at Re = 1685, by means of the 
different forcing wave numbers fcf , are investigated by examining the plotted relative differences in front side impaction 
efficiency in Fig. |4] The peak in relative difference is clear for all three cases, and is found at St = 0.24. So the effect 
of turbulence is largest in the lower region of the boundary stopping mode. The increase in the front side impaction 
efficiency is as expected stronger for the larger forcing length scales. For St = 0.24 the impaction efficiency is almost 
10 times higher for a turbulent flow with an integral scale at kf = 1.5 than for a corresponding laminar flow. The 
case with fcf = 5 also has a dramatic increase in impaction, compared to the laminar, from St « 0.05 on, while the 
change is much less for the kf = 15 case. It should also be noted that the same effects of the turbulence are found at 
the same Stokes numbers for all forcing scales. 



7 




FIG. 2: Relative differences in impaction efficiency for two and tlrree dimensional simulations with forcing at 

kf = 1.5 and fc/ = 5 and Re = 420. 




FIG. 3: Front side impaction efficiencies at Re=1685, for the laminar reference case and the three turbulent forced 

simulations. 



B. Comparison with predictions for large scale turbulence 

In this subsection numerical results will be compared with results obtained from Eq. Q. The expected impaction 
efficiencies, calculated according to Eq. Q, are plotted in Fig. [s] through Fig. 10 The results from the simulations at 



kf = 1.5 (Figures [5] and |6]) are found to be as predicted by the calculations, in the boundary stopping and the classical 
impaction modes. However, in the boundary interception mode for St ^ 0.15, the results from the simulations are 
seen to deviate from what is predicted, especially in the Re = 420 case, shown in Fig. [5j Here, Eq. ([t]) cannot predict 
the somewhat random jumps in impaction efficiency observed at the lowest Stokes numbers. The assumption that 
the impaction efhciency in the inertial impaction mode is not affected by the turbulence, seems to be reasonable. 

Concluding from the results depicted in Figures [5] and |6j the variances of the effective Stokes numbers explain the 
increased impaction efRciencies in the boundary stopping mode. The curves of calculated values i?(77(Stoff )) (shown 
by dashed-dotted lines in the figures) follow the turbulent simulation curves in both figures. 

With kf = 5 and /cf = 15 at Re = 420 (Figures [t] and |9] respectively), the capture efhciency in the turbulent cases 
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FIG. 4: Front side impaction efficiencies for tlie cases with turbulence at Re=1685, plotted as relative difference to 

the laminar. 
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FIG. 5: Calculated impaction efficiency with fcf — 1.5 turbulence at Re = 420, plotted together with the impaction 
efficiencies of case L4 (Turbulent simulation) and case LI (Laminar simulation). 



does not deviate much from the laminar. Therefore, the corresponding scenarios at Re = 1685 (Figures [8]and[l^ 
respectively) are studied when discussing the validity of Eq. ([t]) for fcf = 5 and fcf = 15. 

As seen in Fig. [8| the calculated expected capture efficiency does not match perfectly with the simulated values for 
fcf = 5. In the lower region of the boundary stopping mode, the expected capture efficiency is found to be above the 
values from the turbulent simulation. In the range 0.2 ^ St ^ 0.25, the calculated values match with the simulation 
values, while for 0.25 ^ St 0.90, the simulation values exceed the calculated values. 

Turning to Fig. |10[ the results from the simulation are seen to deviate considerably from the expected values also 
in the boundary stopping mode, up to St « 0.22; the calculated E{r]{^\,c?i)) is too high. For Stokes numbers up to 
St ss 0.15 both the laminar results and the results calculated from Eq. ([7| are found to be higher than the simulation 
results. 

The unpredictable nature of the results seen only in the boundary interception mode (St < 0.15) for the fcf — 1.5 
seems to become more present also in the boundary stopping mode as fcf is raised. This indicates that the behaviour 
of the turbulent eddies close to the boundary layer of the cylinder has a larger effect on particles of increasing Stokes 
numbers as the forcing scale is decreased. 
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FIG. 6: Calculated impaction efHciency with kf — 1.5 turbulence at Re = 1685, plotted together with the impaction 
efficiencies of case H4 (Turbulent simulation) and case HI (Laminar simulation). 
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FIG. 7: Calculated impaction efficiency with ki — 5 turbulence at Re = 420, plotted together with the impaction 
efficiencies of case L3 (Turbulent simulation) and case LI (Laminar simulation). 



C. Particle Clustering 



The mechanism behind particle clustering, or preferential concentration of particles, can be explained as follows. 
The same explanation applies also for backside impaction, as described in the next section. Since the forcing scale 
l{ = Lbox/fcf is a characteristic length for the turbulence, the characteristic eddy time can be defined as 



If 



L 



box 



(19) 



The product Uimsfcf, values of which are found in table |TTj determines Toddy for the different turbulent cases. 

For Re = 1685, the particles exhibited clustering. The Re = 420 cases are not considered, as the phenomenon is 
less apparent here due to the stronger damping of the turbulent eddies. The eddy Stokes number can be introduced 
as 

Steddy = = 3^St. (20) 
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FIG. 8: Calculated impaction efficiency with k{ — 5 turbulence at Re — 1685, plotted together with the impaction 
efficiencies of case H3 (Turbulent simulation) and case HI (Laminar simulation). For St < 0.1, the capture efhciency 

in the turbulent simulation is zero. 
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FIG. 9: Calculated impaction efficiency with fcf = 15 turbulence at Re — 420, plotted together with the impaction 
efhciencies of case L2 (Turbulent simulation) and case LI (Laminar simulation). 



Here, the definition St = Tp/ri has been used. With the characteristic fluid time T[ = D/Uq and the characteristic 
eddy time Teddy given by Eq. (19), Eq. (20 1 can be expressed as 



St, 



'cddy 



D 



Un L 



-St. 



box 



(21) 



It is known that Tn 



"^cddy 1 



St 



eddy 



1, is needed for particle clustering to take place. This is because if 



Toddy S> Tp, the eddy turn-over time is too slow for the centrifugal 'force' to throw the particle towards the wall. 
When Teddy ^ Tp, which is the case for large particles, the particle will not have time to respond to the fast turbulent 
motions and thus will not obtain the acceleration needed. As it is Wmis^f in the prefactor in front of St in Eq. (21) 



that is varying between the different turbulent cases, this product determines for which Stokes numbers clustering 
will be large, namely for those implying Stoddy ^ 1- For the Stokes numbers considered in this work, this implies that 
particle clustering is largest at kf = 15, while the effect gets smaller for decreasing kf. This can be seen in the last 
column of Table [ll| which shows Steddy The discussion of the role of the vorticity in the next section is related to this 
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FIG. 10: Calculated impaction efficiency with fcf — 15 turbulence at Re = 1685, plotted together with the impaction 
efficiencies of case H2 (Turbulent simulation) and case HI (Laminar simulation). 

explanation, since it is the vorticity of the turbulent eddies that gives rise to Steddy, with the magnitude of vorticity 

TABLE II: Some important parameters for the turbulence simulations. Note that the eddy Stokes number in the 

last column is for particles with Stokes number 0.3. 





fcf 




''^rms 


Steddy(St = 0.3) 






15 


0.17 


2.59 


0.13 


Re = 


420 


5 


0.27 


1.34 


0.07 






1.5 


0.42 


0.64 


0.03 






15 


0.24 


3.57 


0.18 


Re = 


1685 


5 


0.28 


1.40 


0.07 






1.5 


0.29 


0.44 


0.02 



D. Backside Impaction 



To further investigate the effects of turbulent eddies on the particles, the impaction on the backside of the cylinder, 
Vhnck, is plotted in Fig. 12 In the Re = 420 simulations, backside impaction takes place only to a small extent 
and is therefore not shown here. For Re = 1685, more general trends in the backside impaction are seen. In all 
cases, the backside impaction is larger than the front-side impaction for St ^ 0.14, and for the turbulence with the 
smallest integral scale 77back > '7front for St < 0.23. The general behaviour for the Re = 1685 cases is that ryback 
increases with increasing Stokes numbers from its initial minimum until it reaches the maximum which is in the range 
0.1 < St < 0.15, from where it decays. 

For a particle to impact on the backside of the cylinder it must first be captured by an eddy bringing it to the back 
side of the cylinder. This could be either a typical von Karman eddy or a large scale turbulent eddy. Being in the 
wake on the back side of the cylinder a particle will typically impact due to an effect very similar to turbophoresis, 
which is encountered in turbulent channel flows. Which particles are impacting on the backside is, in a similar manner 



as for particle clustering (§IVC), dependent on the eddy Stokes number. 

In the laminar cases. Toddy is given by the dimensions of the cylinder, as backside impaction in the laminar cases is 
caused only by the rotational motion of the von Karman eddies in the wake. These wake vortices also play a role in 
the turbulent ryback, but the differences in the turbulent Tyback can be explained by Eq. (21). 
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FIG. 11: Clustering of particles with St = 0.30, from turbulent cases run at Re = 1685. 
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FIG. 12: Back side impaction efficiency for Re — 1685. 



V. CONCLUSION 



Direct Numerical Simulations were run to study the influence of isotropic turbulence on the impaction of particles 
on a cylinder in a cross flow at two different Reynolds numbers, Re — 420 and Re — 1685. The turbulence was 
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simulated on a three-dimensional domain, with forcing at three different wave numbers, i.e. at varying integral scales. 
After reaching a homogeneous steady state, the turbulence was imposed on the main flow domain, with one-way 
coupling between the fluid and the particles. It was checked that for Re = 420 a two dimensional main flow domain 
gave qualitatively the same results as a three dimensional main flow domain. 

The impaction efficiency of particles on the cylinder front side was seen to be greater in the turbulent cases, compared 
to laminar reference cases in the boundary stopping mode. This was found to be related to the statistical variance 
of Stoff, where the effective Stokes number takes into account the fluctuating particle velocities due to turbulent 
velocity variations. The laminar front side impaction efficiency rapidly increases with increasing Stokes numbers in 
the boundary stopping mode, which yields a larger impaction efficiency when turbulence is introduced. The relation 
to the variance of Stcfr was clearest for the largest integral scale, and was less present for decreasing integral scales. 

As the Reynolds number was increased from Re = 420 to Re = 1685, the impact of the turbulence on the particle 
impaction became more dominant. For the smallest particles, in the boundary interception mode, the turbulence 
also lead to changes in the front side impaction efficiency. However, these changes could not be explained by the 
fluctuations in the effective Stokes number. This indicates that the turbulence causes some not fully understood 
mechanism influencing the motion of the smallest particles in the close neighborhood of the boundary layer of the 
cylinder. 

Particle impaction efficiency on the backside of the cylinder in the Re = 1685 turbulent cases was seen to be 
relatively large. This was particularly true for the turbulence with the smaller integral scale, which, due to its strong 
vorticity, had the largest backside impaction efficiency. 

Since the turbulent eddies in the Re = 1685 cases were well sustained throughout the two-dimensional domain, 
the particles dispersed in the turbulent flow exhibited clustering, or preferential concentration. Particle clustering is 
caused by the same mechanisms as backside impaction. Thus, the magnitude of vorticity, and slow dissipation thereof, 
is of crucial importance for the clustering to come into play. 
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